Propagation of the Angular Spectrum in the Up-conversion 
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We study theoretical and experimentally the propagation of the angular spectrum of the light 
produced in the up-conversion process. The connection between the angular spectrum of the fun- 
damental and the second harmonic is derived and measured. We show that even though they are 
connected, it is not possible to directly transfer images from the fundamental to the second harmonic. 
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I. INTRODUCTION 

Parametric up and down-conversion has been subject 
of investigation over the past few decades jj],f|. The 
special case of collinear up-conversion, known as second 
harmonic generation(SHG) has received a great deal of 
attention, because it has become a reliable method for 
generating laser beams in a vast range of wavelengths 
from the green region up to ultra-violet. In both cases, 
up and down-conversion, the transverse properties of 
the field were initially neglected. The interest in the 
transverse properties of the fields involved in the non- 
linear processes is increasing over the past few years ||. 
In most of the theoretical and experimental investiga- 
tions, the fields involved are represented by Hermitc- 
Gaussian modes. More recently, the interest in optical 
modes possessing orbital angular momentum, described 
by Laguerre-Gaussian modes, has found in SHG, a way of 
obtaining modes with higher orbital angular momentum, 
as it is doubled in this process Q] . 

In parametric up and down-conversion processes, the 
converted fields are connected to the pumping fields 
through phase matching conditions. It has been shown, 
for example, that the entangled field of twin photons from 
down-conversion carries the same angular spectrum as 
the pump pi . It was also shown that in stimulated down- 
conversion, the idler beam may carry the same angular 
spectrum as the pump or the complex conjugated of the 
stimulating beam This transfer of angular spectrum 
implies in transfer of images formed by those beams and 
also in transfer of orbital angular momentum |^,^,^|. 

The propagation of the transverse properties of a light 
field, like orbital angular momentum and many others, 
are defined by the propagation of its angular spectrum. 
In this work, we study the propagation of light fields and 
their coupling in the SHG process, taking into account 
their transverse spatial distributions, or in other words, 
propagating their angular spectra. We demonstrated how 
the angular spectrum of the fundamental and the sec- 
ond harmonic field are connected. It is shown that even 



though they are connected, this connection does not al- 
low direct transfer of images from one field to the other. 
On the other hand, as the transfer function is known, it 
is possible to manipulate the second harmonic angular 
spectrum through preparation of the fundamental. 

II. THEORY 

Let us consider the experimental set-up sketched in 
Fig.[|. In the noncollinear case (Fig|l]a) two different 
modes pump a non-linear crystal for generating a third 
mode, the up-converted one. SHG is the collinear version 
(Fig.[I]b), where just one mode pumps the non-linear crys- 
tal and the degenerate collinear up-conversion occurs. 

In the following, we present a quantum theory for 
the up-conversion, based on the formalism developed by 
Mandel and co-workers |3J|. We assume that depletion 
of the fundamental field is negligible in the parametric 
interaction, the second harmonic field is weak and the 
monochromatic and paraxial approximations are used. 
The calculation of the effects of propagation and trans- 
fer of the angular spectrum in up-conversion does not 
require a quantum theory. The same results would be 
obtained with a classical formalism, however the use of 
the quantum hamiltonian and the deduction of the quan- 
tum state for the up-converted field may be helpfull in 
future developements. 

The quantum state of the up-converted field is ob- 
tained in the same way as it is done in Ref. |J for spon- 
taneous and stimulated down-conversion. The difference 
is that here, signal and idler fields are in coherent states 
while the up-converted field is a multimode Fock state 
with n=l: 

= \vac) +C J dq 2 J dq 3 Vi(q 3 - 92)^2(92) |1; ?3>, 

(1) 

where vi(qs — q 2 ) and v 2 (q 2 ) are the angular spectra 
of the input modes with transverse component of the 
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momentum q% — qi and <72, repectively, and 1 ; qr 3 ) is a 
Fock state with n=l for the up-converted mode with the 
transverse component of the momentum qr 3 . 

The transverse intensity profile of the up-converted 
field in a plane situated at a distance z from the crys- 
tal is given by: 



I{r) = {E^(p,z)E^{p,z)), 
where the electric field operator is written as: 



E {+ \p 1 z) = J dqla(q/)exp 



i[qfp 



2k/' 



Performing the calculation we obtain: 



I(r) oc 



dpi Wi(p/)W 2 (pO exp 



(2) 



(3) 



(4) 



For the special case where modes 1 and 2 are the same 
(collinear case), the above equation is simplified to 



J(r) oc 



J dpi W 2 (p/)exp 



(5) 



From the equations above it is seen that the transverse 
field of the second harmonic is given by the square of the 
fundamental, propagated from the crystal to the obser- 
vation plane. 



where q is the projection of the transverse wavevector q 
in the x direction. The angular spectrum is easily prop- 
agated to the plane just before the lens: 

v~(q,z ) oc sine (q~^j cos (q^) exp ^ _i |^ 2 o^ , (9) 

where the paraxial approximation was taken and k is the 
wavenumber. In the plane just after the lens, the angular 
spectrum is given by: 

v + (q,z ) oc J d£ sine (q^j cos ^q^j x (10) 



x exp 



2k 4 



e.xp j ~ij;qS, 



After propagation to the plane inside the crystal and 
turning it back into a field distribution, we obtain : 

W(x, z + f) oc sine (yfxj cos (j§xj x (11) 



x exp 



i jfi(zo - f)x 2 



From the field in Eq. [ll] above, one can calculate the 
final transverse distribution for both fundamental and 
second harmonic fields, in a plane situated at an arbitrary 
distance z from the crystal. For the fundamental it is just 
a matter of propagating the angular spectrum again and 
for the second harmonic one will make use of Eq.|. 



A. Propagation of a double-slit pattern 

Let us consider a plane wave as input to a double- 
slit, as it is represented in Fig.||. Even though this is 
a two-dimension problem, we are going to calculate the 
transverse field distribution at the crystal plane for one 
dimension, for simplicity. After the slits, the beam is 
focused by a thin lens inside the crystal. The lens is 
necessary because the efficiency of the up-conversion is 
low and focusing improves the efficiency || . The outcome 
of this calculation will be used in Eq.|| for calculating the 
final intensity distribution of the SHG field. At the plane 
of the slits, the field is given by: 

W(x,0) = S(x+±a) + S(x-~,a), (6) 

where 



B. Image plane 

The slits are imaged by the lens in a plane after the 
crystal. The longitudinal position of this plane is given 
by the usual relation 
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(12) 



where o is the distance between the slits(object) and the 
lens and i is the distance between the lens and the image. 
According to our scheme, decribed in Fig.|^, o = zq and 
i = f + z D - 

In order to obtain the field distributions at the image 
plane, we write the propagation integral for fundamental 
and second harmonic, starting from Eq.[y], in terms of 
the real parts of the fields: 



S{x, a) 



1 if |x| < f 
elsewhere, 



(7) 



a is the slit width and d is the distance between them. 
The angular spectrum is calculated: 



v(q, 0) oc J dx W(x, 0) exp(— iqx) 
oc sinc(q|) cos(q|), 



(8) 



W(x, z + f + z D )<x fd£ z + f) x 

xexp{ l fc[( 5 ^-^ + ^)e 2 -^+^]} ! 



for the fundamental and 

T(x, z + f + z D )cxfd£ WUt, zo + f)* 



x exp ^2ik 



l 

2-:r 



W 1 



e - Mx 



for the second harmonic. 



(13) 
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The image is obtained when the quadratic part of the 
propagation phase is zero: 



1 

2^ 



2p 



(14) 



This condition is achieved in both cases, when Eq.|l2| 
is satisfied. For the fundamental, it is easy to see that 
when the quadratic term is zero, the field is given by the 
Fourier transform of Wr(£, z + /): 



W{x, z Q + f + z D ) oc 
« / d£ z + /) exp (ik^i 



(15) 



f f 
,, . a - J — 



the image of the slits is amplified by the factor 2y. For 
the second harmonic however, we have: 



T(x, zo + f + zd) oc 
W|(e,z + /)cxp fifc^e 



(16) 



The solution of Eq.[16] is the self-convolution of the 
Fourier transform of the function Wr(£, zq + /), which 
is the double-slit function defined in Eq.[l5|. It is plotted 
in Fig.^ considering the parameters of our set-up. As 
we can see, the image formed by the fundamental is not 
formed by the second harmonic. 



C. Far field 

The intensity distributions for the fundamental and 
the second harmonic in a plane situated at a distance zd 
from the crystal, which is after the image plane, can be 
obtained from the fundamental field inside the crystal, 
given by Eq.|ll|. The analytical solution for the evolution 
integrals is not straightfoward in this case. We have cho- 
sen to solve thisproblem by replacing the real part of the 
function in Eq.|ll| with a sum of gaussians, as an approx- 
imation to solve the propagation integrals. Although a 
closed form of the solution is obtained, resulting function 
is not instructive for the reader. We present instead, in 
Fig.||, a plot this function, for the parameters of our set- 
up, in order to compare it with the experimental results. 



III. EXPERIMENT 

We have performed measurements of the transverse in- 
tensity distribution for the fundamental and the second 
harmonic in the SHG, in different propagation planes af- 
ter the crystal. 

An infrared diode laser, with wavelength centered 
around 845 nm, 150 mW output power is directed to 
a LHO3 (Lithium Iodate) non-linear cristal. SHG takes 
place, producing a beam with wavelength around 425 



nm. Before the infrared reaches the crystal, it is passed 
through a double-slit diffraction screen with slits width a 
= 0.2mm, separation d = 0.4mm and also through a thin 
lens with / = 10 cm focal length. After the crystal, fun- 
damental and second harmonic beams are separated by 
a prism and they are both detected with single photon 
counting modules, based on cooled avalanche photodi- 
odes. The photon counting is necessary for the second 
harmonic beam, because the signal level is strongly de- 
creased by the passage of the fundamental through the 
double-slit. The pulses coming from the detectors are 
sent to photon counters controlled by a computer which 
performs the data acquisition. 



IV. RESULTS AND DISCUSSION 

In a first set of measurements we have placed the slits 
at 12.3cm before the lens, so that the image of the slits 
were formed in a plane situated at 53.4cm after the lens. 
The crystal is always in the focal plane of the lens, so 
that the detection plane is 43.4cm from its center. The 
detectors are scanned in the vertical direction and the 
intensity profile is registered for both the fundamental 
and the second harmonic signals. Fig.|^ shows these pro- 
files. The fundamental presents an extended image of the 
slits, and the amplification factor is in good agreement 
with Eq.|l5|, since d ~ 0.2 mm, 2j- ~ 4 and the distance 
between slits in the amplified image is di ~ 0.8 mm. The 
second harmonic is also in agreement with Eq.|l6|, which 
is plotted in FigJ|. The measurement procedure with a 
finite detection slit smooths the curve. This result shows 
that images can not be directly transferred from the fun- 
damental to the second harmonic. On the other hand, 
the basic information about the image is contained in 
the second harmonic pattern and it might be recovered 
by deconvolution. 

In order to check the evolution of the fields, starting 
from the crystal as it was calculated, we have measured 
the intensity distributions inside the crystal. These mea- 
surements were performed placing a second lens after the 
crystal, imaging the center of the crystal onto the detec- 
tion plane. The results are displayed in Fig.[| These 
patterns are described by Eq.^, where the intensity of 
the fundamental is given by the square modulus of the 
field and the intensity of the second harmonic is given by 
the square modulus of the square of the field. 

In a second set of measurements, the slits were placed 
2cm before the crystal so that no real image was formed 
for the fundamental. In this case, the far field calcu- 
lation must be used. The results are shown in Fig. 0a 
for the fundamental and Fig.^b for the second harmonic. 
As before, we have also measured the intensity distribu- 
tions inside the crystal. The results are shown in FigJ|. 
These curves should be fitted to the square of the field 
in Eq.|ll]. However, we have used gaussian functions for 
fitting the patterns in FigM The parameters obtained 
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with the nonlinear fittings were used in the calculations 
for generating curves in Fig.|]. Comparing the intensity 
profiles in Fig.[I| and the measured distributions in FigJ^, 
a good agreement is obtained. 

For the fundamental, the far field intensity pattern is 
just the same as inside the crystal (FigJ|a), enlarged by 
a scale factor due to the propagation(FigJ^a). For the 
second harmonic, the pattern outside the crystal is dif- 
ferent. The oscillations present two spatial frequencies. 
As a matter of fact, we have fitted the experimental data 
in Fig.^b to a two-frequency interference pattern: 



I(x, zq + f + z D ) oc sine 4 (j^xj X 



(17) 



//i cos [ 2 j^x 



fl2 COS 



The same function could be used to fit the intensity 
pattern of the second harmonic inside the crystal, Fig.||b. 
The reason why FigQb and Fig^b look different is that 
they have different visibilities [i\ and /i2- During the 
propagation the ratio between these visibilities changes 
as it is shown in the theoretical result of Fig.[|b. 



V. CONCLUSION 

In conclusion, we have derived the transverse inten- 
sity distribution for the light beam produced in the up- 
conversion process, as a function of the input pumping 
fields. We have also calculate the intensity distributions 
for the input fundamental and the output second har- 
monic beam generated in the SHG process, when the 
fundamental is diffracted through a double-slit before the 
non-linear crystal. The distributions are calculated for 
the propagation planes inside the crystal, the image plane 
and the far field. It is demonstrated experimentally, that 
images formed on the fundamental are not transferred to 
the second harmonic. It is also demonstrated that the 
double-slit interference pattern transferred to the second 
harmonic, oscillates with two spatial frequencies. The 
experimental data are in good agreement with the calcu- 
lation. 

Financial support was provided by Brazilian agencies 
CNPq, PRONEX, FAPERJ, FUJB and Institutos do 
Milenio-Informagao Quantica. 



M. I. Kolobov, 

Phys. Rev. A 44, 1986 (1991); 

L. A. Lugiato and A. Gatti, 

Phys. Rev. Lett. 70, 3868(1993); 

Gian-Luca Oppo, M. Brambilla and L. A. Lugiato, 

Phys. Rev. A 49, 2028(1994); 

T. P. Grayson and G. A. Barbosa, Phys. Rev. A 49, 2948 
(1994); 

P. H. Souto Ribeiro, S. Padua, J. C. Machado da Silva, 
and G. A. Barbosa, 

Phys. Rev. A 49, 4176 (1994); A. Joobeur, B. E. A. Saleh, 

and M. C. Teich, 

Phys. Rev. A 50, 3349 (1994); 

K. Dholakia, N. B. Simpson , M. J. Padgett and L. Allen, 
Phys. Rev. A 54 R3742 (1996); J. Courtial, K. Dholakia, 
L. Allen and M. J. Padgett, 

Phys. Rev. A 56 4193 (1997); M. Padgett and L. Allen, 
Contemporary Physics 5 275 (2000). 

C. H. Monken, P. H. Souto Ribeiro and S. Padua, 
Phys. Rev. A 57 3123 (1998). 

P.H. Souto Ribeiro, C. H. Monken and S. Padua; 

Phys. Rev. A 60 , 5074 (1999); P.H. Souto Ribeiro, D. P. 

Caetano, M. P. Almeida, J. A. Huguenin, B. Coutinho dos 

Santos and A. Z. Khoury; 

Phys. Rev. Lett. 87, 133602 (2001). 

D. P. Caetano, M. P. Almeida, P.H. Souto Ribeiro, J. 
Hugueni n, B. Coutinho dos Santos and A. Z. Khoury; 
e-print 



A. 



quant-ph/0108101 



A. Mair, A. Vaziri, G. Weihs and A. Zeilinger; 

Nature 412, 313 (2001). 

G. D. Boyd and D. A. Kleinman; 

J. of Appl. Phys. 8, 3597 (1968). 



FIG. 1. Sketch of the experiment. a) Non-collinear 
up-conversion. b) Second harmonic generation. 



FIG. 2. Propagation through a double-slit and a thin lens. 



FIG. 3. Self-convolution for a two-slits distribution. 



FIG. 4. Theoretical intensity patterns, a) Fundamental, 
far field, z = zo + / + Zd- b) Second harmonic, far field, 
z = zo+f + z D - 
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FIG. 5. Experimental intensity patterns, a) Fundamental, 
image plane, z = zq + f + zd- b) Second harmonic, image 
plane, z = zo + f + zd- 

FIG. 6. Experimental intensity patterns, a) Fundamental, 
crystal plane, z = zo + f. b) Second harmonic, crystal plane, 
z = zo + f. 
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FIG. 7. Experimental intensity patterns, a) Fundamental, 

far field, z = zo + f + zn- b) Second harmonic, far field, 

2 = z + f + ZD- 
FIG. 8. Experimental intensity patterns, a) Fundamental, 

crystal plane, z = zo + f ■ b) Second harmonic, crystal plane, 

z = z + f. 
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